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t— ( , Abstract 

Establishing an exact relation for the derivative we show that the eigenvalue flows 
of the hermitean Wilson-Dirac operator obey a differential equation. We obtain a 



On 
O 



ON ■ complete overview of the characteristic features of its solutions. The underlying 

mathematical aspects are fully clarified. 

> ■ 

^ ■ 1. Introduction 

> : 

^ ■ The hermitean Wilson-Dirac operator if is a fundamental quantity in the overlap 

formalism |TI], It has been shown there that the difference of the numbers of its positive 
and negative eigenvalues is related to the index of the massless Dirac operator. This gets 
particularly explicit with the Neuberger operator || . 

The consideration of eigenvalue flows of H introduced in || has initiated a number of 
numerical works on such flows, including studies of the index theorem of topological 
susceptibility 0, of instanton effects 0, and of the spectrum gap with the mass param- 
eter 0. In these works the derivative of the flows has been used f|, [J, considering the 
respective relation as a result of first-order perturbation theory ||. In theoretical consid- 
erations of these flows, as in Section 8 of Ref. 0, details of the behavior at crossing have 
been of interest. Further, questions concerning the eigenvalue flows have also been raised 
in investigations M, |9j of the vicinities of the parameter values —m/r = 0, 2, 4, 6, 8. 

Generally one should be aware of the fact that the considerations of eigenvalue flows of 
H actually imply certain smoothness conditions to hold, which goes beyond the solving of 
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the eigenvalue equation at individual points. Locally this requires adequate properties of 
derivatives and globally that integration provides appropriate solutions. In the case of the 
specific hermitean operator considered, in unitary space fortunately there are theorems 
[ITJ which we can use to settle the first point. To clarify the second point we have to 



develop an appropriate procedure of integration. 

In the present paper we derive a differential equation for the eigenvalues of the hermitean 
Wilson-Dirac operator H and give a complete specification of its admissible solutions. We 
are able to do this in a mathematically well defined way. 

Our developments appear important for a number of problems. In particular, there 
is the ambiguity in the choice of the mass parameter m which affects the counting of 
crossings of flows on a finite lattice. So far, from an upper bound on the gauge field, 
a bounding function has been derived || which e.g. allows to disentangle physical and 
doubler regions. One can hope that the differential-equation properties found here allow 
to extract more detailed information within this respect. In a study of the locality [|TT[] 
of the Neuberger operator [Q] a lower bound for H 2 , also relying on the above gauge-field 
bound, has been important. In such investigations again the differential equation may 
help to get sharper results. In present Monte-Carlo simulations with massless quarks 
(with overlap as well as with domain wall fermions) there are severe problems due to too 



small values occurring for H [fL2 |. The proposals to deal with this include projecting out 



some subspace of small eigenvalues of H ||13|| , constructing forms of H which are better 



behaved around zero [14| , and looking for more suitably gauge-field actions |12] . There are 
now hopes that the differential equation, allowing a more detailed insight, may guide to a 
way out of these difficulties. Studying flows numerically one generally has to interpolate 
between a finite number of points. It appears possible to develop more efficient methods 
for this which use the general properties obtained for the differential-equation solutions. 

In the following we first derive the differential equation for the eigenvalue flows of H 
(Section 2). We then discuss the mathematical properties involved and restrictions due to 
the eigenequation (Section 3). Next we integrate the differential equation and give a com- 
plete overview of its admissible solutions (Section 4). Finally we collect some conclusions 
(Section 5). 

2. Derivation of differential equation 

The Wilson-Dirac operator X/a is given by 

X = ^E V+V, + m + W 7 ,(V, - Vt ) (2.1) 

where (V M ) n ' n = 5 n i n — U^ n 5 n ' tn+ /x and < r < 1. Its property = 75X75 implies that 

H = 75X (2.2) 
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is hermitean. The operator H has the eigenequation 



Hfa = ai<fn (2.3) 

where ati is real and the <pi form a complete orthonormal set in unitary space, as one has 
on a finite lattice. 

Multiplying (|2.3|) by 0)75 one gets <p[^H<j)i = aifyjsfa and summing this and its hermi- 
tian conjugate one has ^{75, H}(pi = 2ai<p\j5<pi. From this by inserting ( |2.2|) with (|2.1|) 
one obtains 

<xi4>h5<i>i = ™- + 9i(m) ( 2 - 4 ) 

where 

9i(m) = ~J2\\VM 2 - ( 2 -5) 
For gi(m) using ||V M 0/|| < ||(V M - t)4>i\\ + \\<pi\\ = 2 one gets 

< gi(m) < 8r . (2.6) 
Further, abbreviating (d«;)/(dm) by d;, we obtain 



<.t. 



^i^z + + <PlH<Pi = $1^1 + c^^p^ (2.7) 

dm dm 

which means that we have 

d, = 0? 7 50i • (2.8) 
Combining (|2.4j) and (2.8) we get the differential equation 

ai(m)ai(m) = m + gi(m) (2.9) 

for the eigenvalue flows of the hermitean Wilson-Dirac operator H. 

3. Requirements for solutions 

In Section 2 actually only continuity of <f>i{m) would have been needed, which can 



be seen repeating the calculations of (|2.7| ) in terms of finite differences. In Section 4, 
analyzing properties of solutions, we shall need gi(m) (at least at certain points) which 
by (|2.5|) implies also the existence of <pi(m). All of this is, however, no problem because 



in the case considered we have derivatives of 4>i(m) up to any order. This follows because 
for our hermitean operator of form H(m) = H(0) + 77275 in unitary space theorems [1C 
apply by which <fii{m) gets holomorphic on the real axis. 

Of the (continuously) infinite number of solutions of ( |2.9|) , specified by integration con- 
stants, only a discrete finite subset occurs for a given H, the selection depending on H. 
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The number of solutions in this subset, being the number of eigenvectors of H, is simply 
the dimension of the unitary space. 

Because of (|2.8|) by the continuity of 4>i{m) also «;(m) must be continuous for all m 
in order that the respective solution a>i{m) of ( |2.9| ) belongs to the admissible subset. 
Since the eigenvectors 4>i(m) have derivatives up to any order this must also hold for the 
admissible ai(m). Because of the continuity required for admissible ai(m) only solutions 
of the differential equation ( |2.9| ) are admitted which exist for all m. 



4. Solutions of differential equation 

Instead of ( |2.9| ) we first consider the differential equation 

/j,(m) = 2(m + 0,(m)) (4.1) 
which by inserting (3i(m) = of (m) becomes ( |2.9| ). Integration of ( |4.1|) readily gives 

rm 

AM = AM) + 2 / dm'M + 9l(™')) (4-2) 



in which particular solutions are determined by the choices of m b and (3i{m h ). These 



choices are restricted here by the fact that one actually wants real solutions of (|2.9| ) which 
meet the requirements discussed in Section 3. 
To get an overview of the properties of (|4.2|) we note that 



dm'(m' + gi(m')) , (4.3) 
where fn is an arbitrarily fixed value, has a minimum at m = vn y if 

m y + gi(m y ) = and gi(m y ) > -1 . (4.4) 

Because of 

H — > 77175 for \m\ — > oo (4-5) 

one has 4>i(m) — > x± with 7sx± = ±X± in this limit. Then in (|2.4|) one gets ai 0^750/ — > 
(±m)(±l) = m so that one obtains gi(m) — >• for |m| — > oo. From this and gi(m) > 
it follows that there is at least one solution m y < of (O). In general several ones with 
m ys < . . . < m yi < m yo < may occur. If there is only one solution of ( |4.4| ) we choose 
nib = iriy. If there are several ones we put nib equal to the m yv related to the lowest 
minimum of ( ^4.3| ) (in case of several degenerate ones picking arbitrarily one of them). In 
this way we achieve that 

rm 

I dm'(m + giiw!)) > for all m . (4-6) 

Jm b 
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In order to get a real solution ai(m) of Q2.9Q which exists for all m, according to (|4.6| ) 



(which takes the value for m = mj,) one has to choose Pi{m^) > in (4.2). At the 



points m yv with minima of (|4.2j ) where (3i(m yv ) > one then immediately sees that 
for the solutions of (|2.9|) one has a minimum of aAm) = + J pAm) and a maximum of 



ai(m) = —y(3i(m). The points where Pi(m yv ) = 0, i.e. the ones with the lowest minima 
of Q4.2D , however, need special consideration. There, to clarify the details at crossing, one 
has (1) to check whether the derivative ai(m) is finite as is necessary in view of ( [2.8|) and 
(2) to disentangle the solutions related to different signs of the square root properly. 
To check under which conditions the derivative a(m) remains finite we note that from 
JD one gets 



• 2 (m + gi (m)) 2 

Oi, = "77 — r (4.7) 

showing that in case of Pi (m) = for some m one must also have m + gi (m) = in order 
that the derivative remains finite. With these relations holding at m one obtains 

af(m) — > 1 + gi(fh) for m — > fn . (4.8) 



Because at the points with Pi(m yi/ ) = envisaged above ( |4.4|) is satisfied, by 
the existence of the derivative af(m yv ) at these points with ai(m y J) = is guaranteed. 
Further, since in ( |4.4j ) one has gi(m yv ) > — 1 it follows that otf{m yv ) > there. Thus, 
because ai(m) should be continuous (as pointed out in Section 3), there must be a crossing 
point of two solutions c^(m), i.e. the solutions ^fJPi(m) from below must continue as 



±w Pi(m) above the zero and have the derivatives 



« i K) = ±/i + 9i K)^ ( 4 - 9 ) 

at the crossing point. 

The asymptotic behavior Pi(m) —>■ m 2 for \m\ — > oo, which implies a(m) — > ±m, should 



be obvious from ( |4.5|) . If desired it can readily be worked out in more detail estimating 



the integral in the form 

P^m) = pi{m h ) +m 2 - m\ + 2 / dm'ft(m') (4.10) 

of ( |4.2|) , which already works using the bound fl2.6|) . 



5. Conclusions 

Establishing an exact relation for the derivatives of the eigenvalues of the hermitean 
Wilson-Dirac operator H we have derived a differential equation for the eigenvalue flows 
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of H. Referring to appropriate theorems also the mathematical aspects have been fully 
clarified. Unambiguous prescriptions for the selection of admissible solutions have been 
given. Integrating the differential equation and analyzing the features of its solutions 
a complete overview has been obtained. Our results appear advantageous for future 
theoretical developments as well as for applications in numerical works. 
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